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The asymptotic behavior of switching time as a function of current for a uniaxial macrospin under 
the effects of both spin-torque and thermal noise is explored analytically by focusing on its diffusive 
energy space dynamics. The scaling dependence (/ — > 0, (r) oc exp(— ^(1 — J)^)) is shown to confirm 
recent literature results. The analysis shows the mean switching time to be functionally independent 
of the angle between the spin current and magnet's uniaxial axes. These results have important 
implications for modeling the energetics of thermally assisted magnetization reversal of spin transfer 
magnetic random access memory bit cells. 
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A spin-polarized current passing through a small mag- 
netic conductor is known to deposit its spin-angular mo- 
mentum into the magnetic system [l] [2]. This in turn 
causes the magnetic moment to precess and in some cases 
even switch direction. This has led to sweeping advances 
in the field of spintronics through the development and 
study of spin-valves and magnetic tunnel junctions (see, 
for example, [3]). A thorough understanding of these phe- 
nomena, however, requires taking into account the effect 
of thermal fluctuations inducing diffusion in the moment 
of a magnetic material. This is of particular experimen- 
tal relevance since spin-transfer effects on nanomagnets 
are often conducted at low currents, where noise is ex- 
pected to dominate. Recent debate in the literature over 
the proper exponential scaling behavior between mean 
switching time and current elucidates how the thermally 
assisted properties of even the simplest magnetic setups 
leave much to be understood [IHS] . 

The standard theoretical approach used to study the 
dynamics of small magnetic structures has been to treat 
the magnetic system as a macrospin in the spirit of 
Brown 10], thus modeling the dynamics as those of a 
magnetic moment undergoing damped precession under 
the action of an effective field. After Slonczewski d , one 
assumes that the magnetic body absorbs spin-angular 
momentum perpendicularly to its orientation. As such, 
spin-torque effects are generally taken into account phe- 
nomenologically by modifying the macrospin's Landau- 
Lifshitz-Gilbcrt (LLG) dynamical equation with the ad- 
dition of the spin-torque term. The LLG equation for a 
normalized magnetization m then reads: 

m = — 7'm X Hcff — aj'm x (m x Hcff ) 

— 7'jm X (m X rip) + 7'ajm x n^, (1) 



'Electronic address: |damele.pinna@nyu.edu| 



where 7' = 7/(1 + a^), 7 is the gyromagnetic ratio, a 
the Landau damping coefficient and j = {h/2e)riJ is the 
spin-angular momentum deposited per unit time with 
r] — (J^ — Ji)/{J-f + Ji), the spin-polarization factor of 
incident current J. The last two terms describe a vector 
torque generated by current polarized in the direction rip. 
The effective field Hcff is immediately obtained once the 
energy landscape of the magnetic sample is specified. In 
the simplified case of a uniaxial monodomain magnet, it 
is characterized by the projection of the magnetization 
onto the uniaxial anistropy axis: 

C/(rrr) = -K{in ■ uk)^ - rrr • H, (2) 

with K = (1/2)MsHkV the uniaxial anisotropy of a 
magnet of volume V and anisotropy field Hk, saturation 
magnetization Ms, applied field H, and the unit vec- 
tor representing the orientation of the uniaxial anisotropy 
(easy) axis. 

Of special interest is the interplay between spin-torque 
and thermal noise effects. Starting from (1), thermal fluc- 
tuations are modeled by considering the effect of a ran- 
dom applied magnetic field with mean zero and variance 
enforced by the fluctuation dissipation theorem. Doing 
this leads to a stochastic LLG equation with multiplica- 
tive noise [TUl [TT] . Experimentally, one is generally inter- 
ested in the magnetization reversal properties of a thin 
fllm magnet where the anisotropy axis is out of plane. 
For currents less than a threshold current (to be defined 
below), the switching in such magnetic systems appears 
to be dominated by thermal-activation related reversal 
processes. Fitting experimental data to theory is com- 
plicated by several factors. First, the validity of the 
macrospin model is questionable for samples typically 
studied (see, for example, [12]), which are so far larger 
than the exchange length of transition metals magnets 
(^ 5 nm). Second, even within the macrospin approxi- 
mation it is not clear what the relevant energetics are for 
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thermally assisted spin transfer torque reversal. This is 
a consequence of the inherently non-conservative nature 
of the added spin-torque term which does not allow the 
construction of a generalized energy landscape suitable 
for Kramer's theory analysis. 

In [S], the authors consider a simplified uniaxial 
macrospin model where the anisotropy and spin polar- 
ization axes were chosen to be collinear. The stochastic 
evolution equation determining the switching properties 
of the macrospin could then be described by a single ID 
uncoupled equation. In this paper we present a parallel 
analysis which we believe to be very useful toward the 
characterization of switching properties in more complex 
macrospin models where symmetry arguments cannot be 
found to simplify the general coupled structure of the 
magnetization evolution equations. 

To elucidate the most basic properties of thermally- 
assisted magnetization reversal under spin-transfer 
torque, we choose to consider a magnetic system solely 
characterized by its uniaxial anisotropy and focus on its 
thermally assisted dynamics in the absence of any applied 
fields. The energy landscape is given by the first term 
on the right hand side of (2). The equilibrium states in 
this landscape are then given by parallel and anti-parallel 
configurations of the magnetization along the uniaxial 
anisotropy axis. We consider our magnetic sample to be 
initially magnetized in thermal equilibrium in one of its 
two possible stable states (i.e. m = +hx or to = —hx)- 
Upon applying a (positive) current, spin-torque effects 
will attempt to drive the magnetization toward align- 
ing with the polarization axis ftp. If the current is then 
switched off again, the magnetization will relax back to- 
ward the easy axis. To capture this behavior, it is very 
helpful to write down the dynamical equation for the 
projection q = m ■ uk of the magnetization along the 
easy axis. For temperatures such that f = K/kT ^ 1, 
the sign of q will specify to a high degree of accuracy to 
which state the magnetization will relax after the current 
is turned off. These "projectional dynamics" read: 

q = a[{nj + q){l- q^) + n.Jqs\ 

+ a'ln,my + ^^{l-q^)oW, (3) 

where / — j /{aHx) is the normalized current and a nat- 
ural time unit t — Hxt has been introduced. If we 
choose ftp as our z-axis, and ft^ oriented somewhere in 
the x-z plane, then s is the projection of the dynamics 
on the axis orthogonal to both the y-axis and uk- To 
simplify the notation, then, and Uz are the compo- 
nents of nx along the original x and z axes. The last 
term appearing in the expression models the effects of 
thermal noise whose stochastic contribution is intended 
according to Stratonovich calculus (indicated by the sym- 
bol o) and is a standard mean 0, variance 1 Wiener 
process PH [T5] . 

In the simplified collinear case, where polarizer and 



uniaxial axes are aligned (riz = 1,71^ = 0), the projec- 
tional dynamics are the same as the regular dynamics 
for the TTiz component of the magnetization, and are de- 
coupled from both rrix and niy . The problem is effectively 
reduced to that of a ID stochastic differential equation 
characterizing the switching behavior over all possible 
currents 

As already hinted, the main complication in studying 
a macrospin's dynamical behavior lies in the generally 
complex structure of its coupled set of stochastic differ- 
ential equations. Equation (3) is precisely a realization 
of such: in the case of collinearity between uniaxial and 
polarizer axes, the evolution equation for = q doesn't 
depend on any of the other two magnetization compo- 
nents. The moment any such symmetry is broken, or 
more complex energy landscapes are considered (such as 
considering biaxial anisotropy scenarios), the full set of 
coupled equations (see (3) in case of general, non-zero 
tilt) must be considered. 

The evolution of a nanomagnet is expected to operate 
in a regime where the nonconservative dynamics act on 
timescales that are much longer than that of the energy- 
conserving part of the hamiltonian. This can be estab- 
lished by noting first that currents / governing the crit- 
icahty of (3) are 0(1) and, second, all terms with the 
exception of the first in (1) are of order a or higher. As 
such, for a ^ 1, the two evolution timescales are sep- 
arated. Trajectories then should remain close to stable 
periodic Hamiltonian orbits, and slowly drift from one or- 
bit to another due to damping, STT and thermal noise. 
In [5] , the authors addressed the problem of broken sym- 
metries due to the presence of tilt by averaging equation 
(3) over constant energy trajectories. In doing so, they 
again retained a ID stochastic differential equation and 
solved its associated mean first passage time problem. 
This approach, though, is just a particular example of a 
more general analysis which can be applied to macrospin 
dynamics. 

Take as an example the form of the energy landscape in 
(2) in the absence of applied magnetic fields and ask how 
it evolves in time. Taking a time derivative and using the 
definition of q, one finds: 



Combining (3) and (4) and averaging the dynamics over 
constant energy trajectories leads to a diffusion equa- 
tion in energy space [U [T31 |T3]. Such an equation is 
one dimensional and can be defined independently of the 
complexity of the energy landscape for any generalized 
macrospin system. In our particular case, the averag- 
ing is very straightforward to accomplish since constant 
energy trajectories (a — > 0,/ — > 0,^ — > oo) are simple 
circular librations around the uniaxial anisotropy axis 
with q ~ sgn((7)-\/|ej. Proceeding with the averaging and 
noting that (g^s) = {qrriy) = due to the geometrical 
structure of the constant energy orbit, we find for states 
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starting antiparallel to the uniaxial anisotropy axis: 



|6| = 2a^\e\{l - \c\)W\e\ - In^) + 2^-|6|(l - |e|) o W. 

(5) 

Two fundamental properties of this energy space evo- 
lution are immediately seen. The first is that In^ = 1 
behaves as a crossover value above which all dynamics 
prefer evolving towards the parallel state independently 
of noise. Accordingly, this value of /n^ is hereafter de- 
noted as the "critical current". Analogously, the ther- 
mal regime is characterized by sub-critical current values 
{Iriz ^ 1). The second is that the angular tilt between 
uniaxial and polarizer axes factors in only through the 
term multiplying the current /. As such, as the angu- 
lar tilt is varied, dynamics are going to differ only up to 
a rescaling of the current / and, beyond that, all other 
properties of the dynamics will be functionally identical. 
This was one of the main results shown in [9^ where the 
authors averaged (3) along constant energy trajectories. 

The thermally assisted properties are then easily un- 
derstood by studying the mean first passage time over 
the effective barrier to the dynamics of |e| analogously 
to what has been done in the literature. Instead, we 
consider another method (initially due to Friedlin and 
Wentzell [H]) to extract the exponential scaling behav- 
ior of the mean first passage time without resorting to the 
full mean first passage time problem. Formally, given a 
Langevin equation x — f{x) + g{x) o W, the Wiener pro- 
cess W can be solved at each instant in time. At each 
such instant, the probability of a Gaussian-distributed 
random kick can be written as: 



P[W]\t=exp 



cxp 



(6) 

The probability of specific realization of a stochastic path 
will then simply the product of the probabities of all the 
random kicks that the path witnessed. As such the prob- 
ability of a path starting in some point A and ending in 
B in time T will be: 



P[A B]^ exp 



2 V gix) 



(7) 



The most likely trajectory starting in A and ending in B 
in time T will be that which maximizes its probability of 
realization, or conversely, that which minimizes the F-W 
action: 



'S'fw — / dt- f - — J-Y^ 
Jo 2 V .9(2;) 



(8) 



Finding the action of such a trajectory is then a standard 
Hamilton- Jacobi problem: 



A5fw — Sfw{B) — Sfw{A) 



dx. (9) 



We are interested in the F-W action of the most prob- 
able trajectory starting in |e| = 1 and ending at the sad- 
dlepoint of the deterministic dynamics: |e| = (we spe- 
cialize to the case = 1 since all other cases are imme- 
diately reobtained be rescaling). Plugging (5) into (9), 
one finds that: 



A^FW — C 



/2 



(10) 



The exponential scaling between mean switching time 
and current is then r oc exp(^(l — 7)^) 

To verify this result, we numerically solved the full 
stochastic LLG equation (1) by using a standard sec- 
ond order Heun scheme to ensure proper convergence to 
the Stratonovich calculus [17] . As mentioned above, the 
nanomagnet is initially in thermal equilibrium at a tem- 
perature T in one of its two possible stable states at the 
start of the simulation (t = 0), when the current is ap- 
plied. At each time step, the random kicks are generated 
with strength given by 
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a 



2if(l + a2) 2(1 + a2)^ 



(11) 



which is obtained by imposing the fluctuation-dissipation 
theorem in the absence of applied currents. Such a diffu- 
sion constant is valid as long as the magnetic sample is 
in equilibrium with the thermal bath and effects due to 
Joule heating are neglected. 

As already discussed, upon breaking the symmetry of 
the uniaxial problem by introducing a tilt between uni- 
axial and polarizer axes, the scalings should remain iden- 
tical up to a substitution / — > /n^. Furthermore, analo- 
gously to the coUinear limit, one can again define a criti- 
cal current 1^ as the current above which all initial mag- 
netic configurations will switch deterministically. This 
has already been derived elsewhere in the literature [7] 
and can be shown to give /c = 1/ cos(a;), where oj is the 
angular tilt between uniaxial and polarizer axes. This 
in turn automatically implies that Uzl = I / Ic- In other 
words, the thermally activated dynamics are the same 
for all angular separations up to a rescaling by their re- 
spective critical current. This result is shown in Fig. [T] 
It is of interest to point out how the effect of angular 
tilt in thermally activated spin-transfer reversal is quite 
different from what one gets in studying the similar phe- 
nomenon of field driven reversal. It has been shown, in 
fact, that upon introducing an angular separation be- 
tween easy and applied field axes, the scaling exponent 
changes abruptly in value from 2 for the coUinear case to 
3/2 [IH]. 
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FIG. 1: Comparison of the mean switching time versus cur- 
rent in the sub-critical low current regime for different val- 
ues of angular separation between polarizer and easy axes, 
oj, ranging from to 7r/4. All critical currents have been 
rescaled to 1. Times are shown in units of (s • T) where T 
stands for Tesla: real time is obtained upon division by Hk- 
Black curve is theoretical fit from [J where the exponential 
scaling behavior r oc exp(^(l — J)^) was obtained by other 
means. It is apparent that, up to a current rescaling, the 
mean switching time vs applied current is independent of cu. 



These results have important implications for the anal- 
ysis of experimental data in which measurements of the 



switching time versus current pulse amplitude are used 
to determine the energy barrier to magnetization rever- 
sal. Clearly, the correct asymptotic scaling form is es- 
sential to properly determine the energy barrier to re- 
versal. The energy barrier, in turn, is very important 
in assessing the thermal stability of magnetic states of 
thin film elements that are being developed for long term 
data storage in STT-MRAM. Further work should ad- 
dress how these results extend to systems with easy plane 
anisotropy and situations in which the nanomagnet has 
internal degrees of freedom, leading to a break down of 
the macrospin approximation. The method exposed al- 
lows to treat analytically even more complex dynamics 
where other anisotropy contributions wish to be tken into 
account. To keep the exposition as linear as possible, we 
have refrained from solving the full mean first passage 
time problem associated to equation (5). Doing so is 
straightforward and proceeds exactly along the lines of 
what was shown in [5]. The one advantage of solving 
the full first passage time is that one is capable of ex- 
tracting information on the exponential prefactor for the 
reversal time scalings. We will address the application of 
this methodology to the full biaxial anisotropy model in 
the presence of spin-torque and thermal noise in a future 
publication. 
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